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Abstract: In an observed semi-Markov regime, estimation of transition rate of regime switching 
leads towards calculation of locally risk minimizing option price. Despite the uniform convergence of 
estimated step function of transition rate, to meet the existence of classical solution of the modified 
price equation, the estimator is approximated in the class of smooth functions and furthermore, 
the convergence is established. Later, the existence of the solution of the modified price equation is 
verified and the point-wise convergence of such approximation of option price is proved to answer 
the tractability of its application in Finance. To demonstrate the consistency in result a numerical 
experiment has been reported. 
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1 Introduction 


Among the researchers of financial market modelling and derivative pricing, regime switching econ¬ 
omy is a popular choice, owing to its simplicity and mathematical tractability, along with its ability 
to incorporate fluctuations of market parameters. See DiMasi et al 1994 ; Mamon &: Rodrigo 
[2005 ; Joberts & Rogers 2006 ; Deshpande & Ghosh |2008| ; Basak et al 2011| ; Ghosh & Goswami 


[ 2009 and references therein for more details. Broadly speaking, regime switching economy refers 
to a class of mathematical models of financial securities where the market parameters such as ex¬ 
pected growth rate, volatility, interest rate etc are assumed to evolve as a stochastic process with 
finite state space. Therefore, for obvious reason, a fair price of a contingent claim under such 
market assumption, depends on the law of this finite state process. In practice, the law of such 
regime switching is not known a priori. If the regimes are observed, then one can perhaps estimate 
the transition rate of the regime switching and calculate an option price using the estimator. This 
approach would be satisfactory provided the calculated price is ‘close’ to the theoretical price in 
some sense. In this paper we examine validity of such approach. Since the market, under discussion, 
is incomplete, there are multiple fair prices. We consider the unique locally risk minimizing price 

and often do not mention explicitly. 


as m 


Schweizer 2001 


In this paper we consider semi-Markov regimes which subsumes the popular Markov regimes as 
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well. It is known that unlike Markov case, the transition rate in this generality cannot be written 
as a constant matrix. The rate of semi-Markov process turns out to be a matrix valued measurable 
function on [0,oo). Nevertheless, it is possible to construct a sequence of MLE of certain finite 
dimensional projections of the transition rate function and then establish the almost sure conver¬ 


gence of the sequence to the transition rate function Sansom Sz Thomson 

[] [I^ 


2001 ; Ouhbi Sz Limnios 


On the other hand, it appears that the price function of a European option satisfies a 

for more details) in which the transition rate 


differential equation (see Ghosh &: Goswami 2009 


appears as a parameter. It is shown in Goswami et al ] (M5l that this equation admits a classical 
solution if the transition rate satisfies certain conditions including differentiability. Therefore it is 
important to ensure two things. Firstly, if the transition rate function is replaced by a smoothed 
version of a member of the sequence of MLE, in the differential equation, does the equation still 
admit a classical solution? Secondly, if one gets a sequence of approximated price in the above 
manner, does it converge to the theoretical price? We answer both of these in the affirmative. We 
also illustrate this convergence result by one numerical experiment. 


The rest of this paper is arranged as follows. In Section 2 the mathematical model of the financial 
market and the pricing equation is briefly presented. A particular smooth approximation of the 
transition rate function is constructed in Section 3. It also contains the convergence result of the 
smooth approximation. In Section 4, we show that using the sequence of estimator of transition 
rate function, one can construct an approximating sequence of the European call option price which 
converges to the true price. Finally some numerical experiments are reported in Section 5. 


2 Market model and pricing 


Here we consider a market consisting of only two assets among which one is a locally risk-free asset, 
such as a money market account and another is a risky asset, such as stock. Let {W}t>o be a 
semi-Markov process on a finite state space X = {1,2,...,0} with transition probabilities from 
state i to state j as pij and conditional holding time distributions given the present state i and 
next state j as F{-\i). Assume that (H, P) is the underlying complete probability space which 
also contains a Wiener process {Wt}t>o that is independent to {Xt}t>o- We model the hypothetical 
state (assumed to be observable) of the market by {Xt}t>o and assume that the bank interest rate, 
volatility coefficient, and the growth rate evolve as functions of Xt and are denoted as r{Xt), cr{Xt) 
and lJi{Xt) respectively. We denote the consecutive transition times of Xt by {rn}n>o with Tq = 0 
and the holding time at t by 1*. The price of risk-free and risky assets are denoted by {Bt}t>o and 
{St}t>o respectively and are given by 


Bt 

dSt 


exp( / r{Xu)du), 

Jo 

St{p{Xt)dt + a{Xt)dWt), So >0 


( 1 ) 

( 2 ) 


where a is positive valued and r, /r are nonnegative maps. Let J-t be the filtration of P generated 
by Xt and St and satisfying the usual hypothesis. In Ghosh X Goswami 2009 it is shown that 


this market model does admit an equivalent martingale measure, thus it is arbitrage free under 
admissible strategy. It turns out that the above market is incomplete. Therefore, no-arbitrage price 
of a derivative might not be unique. In this paper we consider locally risk minimizing option price 
which exists uniquely Ghosh &: Goswami 2009 . For more details about this pricing approach in 


a more general setup we refer Schweizer 2001| and the references therein. 

In order to obtain a representation of the option price we further assume that 
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(Al) (i) F(y|i) e C2([0,oo)) Vi 
(ii) F{y\i) < 1 Vi and ?/ > 0 
(ui)(pij) is an irreducible probability matrix. 


We denote the instantaneous transition rate function as \ij{y) which is given by 

-P(^T„+1 =j,y< ^T„+1 <y + At/I = i, > y) 


lim 


Ay 


Thus Xijiy) = Pij 

In this paper, as an example, we consider a European call option, on St with strike price K and 


maturity time T. Then the contingent claim {St — K)^ is Ft measurable. It is shown in Ghosh 


& Goswami 2009| and Goswami et al 2015| that the locally risk minimizing price of this claim 


at time t is a function of t,St,Xt and Yt and that, ip say, is the unique solution in the class of 
functions of at most linear growth of the initial boundary value problem 


defined on 


d d d 1 d'^ 

—V7(t, s, i, y) + s, i, y) + r{i)s—<p{t, s, i, y) + -a‘^{i)s‘^-^(p{t, s, i, y) 

'^Xij{y)[Lp{t,s,j,^) - ip{t,s,i,y)] = r{i)p{t,s,i,y), 
j¥=i 


X A X {0,T)\y £ (0,t)}, 


( 3 ) 


( 4 ) 


V := {{t,s,i,y) £ {0,T) x 
with boundary conditions 

lim(^(t, s, T, y) = 0, VtG[0, T], 
s4,o 

^{T,s,i,y) = {s - K)^; s G M+; 0<y<T-, i = l,2,...,9. (5) 

In Goswami et al [2015| it is shown that the above problem is equivalent to a Volterra integral equa¬ 
tion of second kind which can be comfortably solved numerically using a step by step quadrature 
method. 


3 Approximation of transition rate 

We augment the semi-Markov process Xt with the holding time process Yt to obtain a process 
{X,Y) = {Xt,Yt)t>o, which is clearly Markov. 

Because of the deterministic nature of {Xt,Yt) during every interval (T„,r„_|_i), the dynamics of 
{Xt, Yt) can be described by a discrete time Markov process (A„, Yn)n>o where we define Xn '■= Xt„ 
and Yn := Yt^- = Tn — T„-i. Thus the augmented Markov process can uniquely be specified with 
initial distribution as P{Xq = j) ■ = p{j), Ytq = 0 and the following semi-Markov kernel, 

P{Xn+i =j,Yn+i <y\Xo,Xi,...,Xn,Yi,...,Yn) := px„jF{y\Xn) {a.s.) (6) 

for all y G M and 1 < j < 9. We denote Qij{y) := pijF{y\i) Vi / j, Aj(y) := Xij{y), 

^i{y) ■= fo Xi{u)du. We see, F{y\i) = Yfj=iQij{y) and also 

= /(^N) = X] = (1 - ^(''^10) X] 
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By solving the above ODE for E, one obtains 


ln(l — F{y\i)) = —Ai{y), which implies F{y\i) = 1 — exp(—Aj(y)). 


Hence 

e 

^Qij{y) = 1 - exp(-Ai(y)). (7) 

t=i 


Consider a history of augmented Markov process censored at fixed time r, 




where AA is the number of transitions before time r and LA := r — Tjv^ is the backward recur¬ 
rence time. The associated log-likelihood function is maximized to obtain the maximum likelihood 
estimator(MLE) of the transition rate function, The likelihood function for 77(r) is 

e Nr-l 

L{T)=p{Xo){l-Y,Qx^j{Ur)) n Px,x,^J{Yi+i\Xi). 

l=l 1=0 

Thus from Q, 

N^-l 

p{Xq)-^L{t) = exp{-Ax^^{Ur)) exp{-Kxi{Yi+i))^Xi,Xi+^{Yi+i). 

1=0 

Then we consider log-likelihood as 


N-r-l 

1{t) := log{p(A:o)"^L(r)} = ^ (log Ax,,x,+i(Ti+i) - Ax,(Ei+i)) - Ax^^iUr). 

1=0 

We consider {vk)o<k<M-iy a regular subdivision of [0,r] with step A,- = ^ and M = where 

a > 0, to define for i 7^ j G T and y > 0 

M-l 

\j{y) (y)’ (^) 

k=0 

where A^-fc = Xij(vk)- After replacing A by the step function A* in the expression of 1{t) one obtains, 

M-l 

EE los Ajjffc Xijk'Uiii'), (9) 

i,j£S k=0 


where 

Vik ■= ^ (Yl+l A Vk+l - Vk)l{i}x(vk,oo)iXl, Yl+i) + {Ur A Vk+1 - Vk)l{i}x(v^:,oo){^Nr,Ur), 
1=0 


and 


Nr-l 

dijk ■= ^{i}x{j}x{vk,v,,+i]{Yi, Xl_^_l,Yl^l). 

1=0 
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Hence the estimator of Xijk which maximizes the functional in Q, is given by 



dijk/Vik if Vik ^ 0) 
0 otherwise. 


Thus we obtain an estimator of X*j{y), given by, 

M-l 

^ijk^{vk,Vk+i](y) + ^b'ol{0}(y)- 

k=0 


( 10 ) 


( 11 ) 


We have the following result. 

Lemma 3.1. Fix a G (0,1/2). Under the assumptions of (Al), the estimator r), is uniformly 
strongly consistent for Xij{■), on [0, T] in the sense that 

max sup |Aij(y, r) — Ajj(y)| —>• 0, almost surely, as r —)• oo. 
ye[o,r] 


Proof. If (Al) holds, one can directly derive from the definition of Xij{y) that Xij{y) = 

Now by summing over j € X \ {i} and then integrating on [0, y) both the sides, one obtains 
Ai{y) = — ln(l — F{y\i)). Thus under assumption (Al) Aj(y) is in C^([0, oo)). Hence, for any i ^ j, 
Xij : [0, oo) —)• [0, oo) is well defined and continuously differentiable. Now the rest of the proof 

[T^. □ 


follows from Theorem 1(b) of Ouhbi & Limnios 


In view of the above Lemma, we fix a G (0,1/2) for the rest of this paper. Note that A(-,r) 
is a step function and thus it is discontinuous. We aim to obtain an approximation of A in the 
class of smooth functions so that the approximation can be used to obtain an approximated price 
function by solving appropriate system of differential equations. To this end we consider B^- 
spline interpolation of A(-,t). We first extend it on M by assigning X{y,T) = A(0,t) Vy < 0 and 
X{y,T) = A(Ti,r) Vy > Ti := arymax^gj-^,-] A(y,r) and setting Vk '■= kA^-, for all A; G Z. Set 


Bliy) 


{Vk+2 - Vk){Vk+l - 

( jy - Vk){vk +2 - y) {vk+3 - y)iy - Vk+i) A , , 

\{Vk+2 - Vk+l){Vk+2 - Vk) {Vk+2 - Vk+l){Vk+3 - Vk+l) J 

I (^fc+3 - y)^ . / ^ 

^(Ufc+S - Vk+2){vk+3 - Vk+l) 


For each i j, we interpolate the data points (r^, Ajj(;._i)) and we denote the spline interpolation 
of Xij as Xij which is given by 


OO 

hj{y,T):= ^ Xij{vk+2,T)Bl{y) = ^ kj{k+i)Bl{y) 


( 12 ) 


k=—oo 


k=—oo 


Since for each i / j, Xijk is non-negative for every k, from the property of B-spline, Xij is a 
non-negative function. 

Lemma 3.2. Under (Al) 

max sup |Aij(y,r) - Aij(y,r)| 0 
ye[o.r] 

as T ^ oo. 
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Proof. From Theorem 6.6.4 of Kincaid &: Cheney 2010 , for each i ^ j, supyg[o,T] “ 


Aij(?/,r)| < 2a;[o,T](Aij(-,r); A.r) where u!iif,d) := sup{|/(f) - /(s)||f,s € I, |i - s| = <5} is the 
modulus of continuity. Again it follows from Lemma |3.1[ that for a given e > 0 there is a A such 
that P{N < oo) = 1 and for t > N, maxj^j supj^g[o,T] |A*j(y,T) — \ij{y)\ < s. Thus for t > N, 
‘^[o,T](Aij(-, t); At) < W[o,r](Aij(-); ^r) + 2e. Since Ajj(-) is continuous, its modulus of continuity 


converges to zero as A^ 
result. 


0. Therefore, we get as r —)> oo, r); A,-) 


0 . 


Hence the 

□ 


Lemma 3.1 and 3.2 lead to the following result. 


Theorem 3.3. Under (Al) 


max sup \\ij{y,T) — Xij[y)\ —>■ 0, a.s., as t ^ oo. 
[0,T] 


4 Approximation of price function 

It is shown in the previous section that a convergent sequence of smooth approximations of transi¬ 
tion rate function can be constructed using a combination of non-parametric MLE and B 2 spline. 
Having this result, it is tempting to solve Q-([^ with the smooth approximation of transition rate 
to obtain an approximation of price function. Needless to mention that such approximation is 
reliable only when certain continuous dependency of the solution on the transition rate function is 
established for the concerned initial boundary value problem. Such result is not readily available 
for the non-local degenerate type of parabolic PDF, we consider here. In this section we establish 
the convergence of such approximate price function to the true price. For the sake of preciseness 
we propose the following definition. 

Definition 4.1. Let X : X x X x [0, 00 ) — )• [0, 00 ) an approximation of the transition rate function 
X and ip the solution of the problem (©-(©. If (I)-©, after the function X is replaced by X, admits 
a unique classical solution ip, then ip is called the TEA (Transition rate based approximation) of (p 
with parameter A. 


Goswami et al 


2015 


presents a set of fairly general sufficient condition on A, which ensures 
existence of a unique classical solution of the modified equation. 

Lemma 4.1. Under (Al), the TBA of p (as in (I)-®; with parameter X{-, t) (as in ([l2|) )exists 
for sufficiently large value of t. 

Proof. We define A(y,r) := A^-(y, r), Ai{y,T) = Xi{y,T)dy, F{y\i,T) := l-exp(-Ai(y, r)), 

and pij := l(0,oo)(Ai(y, T))dF(y\i, r). If we can show that there is a tq > 0 such that for 

hold i.e., (i) A(-,t) is in C^([0,oo)), (ii) 


'0 Ai(y,r) 

T ^ To j the sufficient conditions in 


Goswami et al 


2015 


lim^^^oo A.i(y, r) = 00 for each i, and (iii) pij is irreducible, then the existence of TBA of p follows 
from Theorem 3.2 of Goswami et al ] (m^ . 

We note that from the property of B-spline, (i) holds. Again since for each i, limy_,.oo Aj(y) = 
3 — ln(l — F{y\i)) = 00 , from Theorem 3.3 there is a tq > 0 such that for r > tq, (ii) holds. 


lim. 


We first note that to prove (iii), it is sufficient to show that if pij > 0 for some i j, then pij is 
also positive. Now if pij > 0, Xij{y){= Pijiz^pj^) is not identically zero. Hence using Theorem 
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3.3 


there exists a ri > 0 such that for any r > ri, Xij{y,T) is not identically zero. Thus 

Pij = [ 1(0,oo) (A»(j/, T))dF{y\i, r) 

Jo Aj(y,T) 

1*00 

= / Kj{y,r)e~^J^''^Uy 

Jo 

> 0 . 


□ 

We state the main result below. 

Theorem 4.2. Under (Al) for a large value of t, let, (fU he the TEA of ip with parameter A(-,t). 
Then, ip'^ converges to ip point-wise, as r —>■ oo. 


In order to prove this theorem, we need the the following lemma. 

Lemma 4.3. Let {St}t>o o,s in (§. Then 

E[£ Stdt]<^{e^^-l) (13) 

for some positive constant C. 


Proof. We know, Q has a closed form solution given by 

1 


St = So exp 
We introduce the following constants 


f {pt{X^)- a\Xu)}du+ / a{X^)dWu 

0 ^ Jo 




c := max p,{i) — -cr^(i) /• , d := max{iT^(i)}. 


,2/ -N 


iex 


Clearly, 


We observe, 


St < So exp(ct) exp a{Xu)dW^ 


: °° pTnAt 

a{Xu)dWu = V / a{XT„_MW, 


J Tn—lAt 




(14) 


(15) 


n=l 

Hence, the conditional distribution of a(Xu)dWu given is normal with mean zero and variance 

oo 

At — Tn-1 A t), 

n=l 

where is the filtration of X generated by X = 
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Now using the formula of variance of a log-normal random variable, we get, 


E 


exp / a{Xu)dW^ 


= E 


= E 


< E 


E 


exp / a{Xu)dWi, 




■X 


exp ( - ^ a^{XT^_^){Tn At- r„_i A t) 


n=l 
, oo 


exp - A t - Tn-i A t) 


n=l 


= exp ( -t 


Using above inequality we obtain from (15), 


E{St) < 5oe"*ei* = 


Since St is nonnegative we apply Tonelli’s theorem and the above relation to get, 


E 


Sfdt 


Ljo 


EStdt 


< So [ 

Jo 

- l) . 


(c+i) 


□ 


Proof of Theorem 


4.2 


We define the difference of the functions (/? and its TBA (fE to be. 


s, i, y) := (/?(t, s, i, y) - (fE{t, s, i, y). 


Now, by considering the initial boundary value problems satisfied by ip{t,s,i,y) and cp'^{t, s,i,y) 
respectively, one obtains directly that satisfies the following initial boundary value problem, 

d d 1 

—'ip{t,s,i,y) + r{i)s—'ip{t,s,i,y) + s,i,y) + '^\ij{y){'il){t, s, j,Q) - i/jit, s,i,y)) 

= r{i)'il;{t,s,i,y) {t, s, j,0) - {t, s,i,y)), 

j¥=i 


defined on 


with conditions 


V := {{t,s,i,y) € (0,r) x M+ x T x {0,T)\y G (0,t)}, 


lim?/)(t, s, i, y) = 0, Vte[0, T], 
sto 


s,i,y) = 0, s G M"' 
We rewrite the above system of equations as. 


0 < y < T; z = 1,2,- • • ,6». 


—s, z, y) + Cipit, s, z, y) = s, z, y) - f{t, s, z, y) 


(16) 


















where 


{C'ilj){t,s,i,y) := 


r-\ 9 1 2 /•^ 2 

'■W'ai + 2- M» ^ 


ip{t,s,i,y) + '^Xij{y){i^{t,s,j,0) - ipit, s,i,y)), 
j¥=i 


and 


r{t,s,i,y) ■■= ^(Aij(y) - Xij{y,T)){ip^{t,s,j,0) - s,i,y)). 

j¥=i 


Note that, L is the infinitesimal generator of {St,Xt,Yt) satisfying 

dSt = St{riXt-)dt + a{Xt-)dWt), 

where Xt is a semi-Markov process with transition rate Xij{y) and Yt is the holding time process. 
Then, using Feynman-Kac formula, 


/ exp - / r{Xu)du\rf {t')dt'\St = s,Xt = i,Yt = y], 


(17) 


where. 


From Theorem 


3.3 


y^{t) = rit,St,Xt,Yt). 

we have the uniform convergence of Xij{y,T) to Xij{y) in y G [0,T], almost 
surely. Or, in other words we get, for any given e(> 0), 3N s.t. P{N < oo) = 1 and for t > N, 

\Xij{y,T) - Xij{y)\ < s VyG[0,r]. (18) 

Using Theorem 3.2 of Goswami et al ] there exists constants ki and k 2 such that 0 < 

s,i,y) < ki + k 2 S for all i,y and t> N. Hence from (18), 

r{t,s,i,y) < e{ki + k 2 s). 


Therefore, by applying Lemma 4.3, we assert that for all large r, \r]'^{t,u})\ is dominated by a fixed 
integrable function. Thus using dominated convergence theorem, we can see 

lim'il;'^{t,s,i,y) = lim E[J^ exp^-^ r{Xu)d't?j r]'^{t')dt'\St = s, Xt = i,Yt = y] 

= E[J^ exp^-^ r{Xu)dv^limr]'^{t')dt'\St = s,Xt = i,Yt = y] 

= 0 . 


Hence, (p'^ converges to ip point-wise, as r —>■ oo. □ 

Remark 4.1. In this section, as an example, we consider a European call option to illustrate the 
the convergence of TEA of price. It is important to note that one can have similar convergence 
results for other derivatives such as put option price, barrier option price, compound option price 
etc. as the prices satisfy similar system of PDEs. 
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5 Numerical experiment 

We firstly illustrate the convergence result as in Lemma |3.1| through considering an example of a 
semi-Markov process and estimating its transition rate over different time scales of the history of 
the process. In order to compute and present the actual estimation error, we simulate the example 
having three hypothetical states {1,2,3} and the holding time distribution as 

f{y) ■= for y>0, 

and we choose the transition matrix as 

/ 0.0 0.1 0.9 \ 

(Pij) = 0.4 0.0 0.6 . 

\ 0.7 0.3 0.0 / 

Therefore, the theoretical transition rate function for any y > 0, is given by, 

Aii (y) = Pij . = Pij yi,j ^ S and i / j. 

M-F(y|^) ■'y + l 

The convergence as r —)• oo is illustrated by plotting two different norms of the difference between 
A and A* on [0, T] where T = 4 for many different values of r. Although the convergence result 
is established only for sup norm, but we find it interesting to illustrate the norm also. In the 
Figure r variable is taken on the horizontal axis and the computed norms of error is plotted on 
the vertical axis. In the plots, logarithmic trend lines are fitted to visualize the trends. To illustrate 



Figure 1: Convergence of MLE 

the convergence of TEA of a European call option price function we consider a market where the 
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volatility and instantaneous interest rate at each regime are given as follows 


(r(l),r(2),r(3)) = (0.3,0.6,0.7) 

(a(l),a(2),a(3)) = (0.2, 0.2,0.2). 

The strike price X is 1 and the maturity T is 1 unit. Figure shows convergence of TEA of price 
function through the plot of sup norm of the error for the estimation, computed as 

Error := ||V’^(0, •,i,0)l[o,5]||snp- 

The norm of "(/^^(O, •, i, 0)l[o_5] is also plotted and finally a logarithmic trend line is added to each 
of the data series. 



Figure 2: Convergence of approximation error 


Acknowledgement: The authors are grateful to Jayant Deshpande and Mrinal K. Ghosh for 
very useful discussions. 


References 

Basak G. K., Ghosh Mrinal K. and Goswami A., Risk Minimizing Option Pricing for a Glass of 
Exotic Options in a Markov-Modulated Market, Stoch. Ann. App. 29:2(2011), 259-281. 

Deshpande A. and Ghosh M. K., Risk Minimizing Option Pricing in a Regime Switching Market, 
Stoch. Ann. App. 26(2008), 313-324. 

DiMasi G. B., Kabanov Y. and Runggaldier W. J., Mean-Variance hedging of options on stocks 
with Markov volatility. Theory Probab. AppL, Vol. 39 (1994), 173-181. 


11 





Ghosh M. K. and Goswami A., Risk Minimizing Option Pricing in a Semi-Markov Modulated 
Market, SIAM J. Control Optim. 48(2009), 1519-1541. 

Goswami, A., Patel, J., Shevgaonkar, P., A system of degenerate non-local parabolic PDE and 
application, 2015. arXiv:1506.01467 [math.AP]. 

Joberts A. and Rogers L. C. G., Option pricing with Markov-modulated dynamics, SIAM J. Control 
Optim. 44(2006), 2063-2078. 

Kincaid D. and Cheney W., Numerical Analysis: Mathematics of Scientihc Computing, American 
Mathematical Society; 3rd Revised edition (2002). 

Mamon R. S. and Rodrigo M. R., Explicit solutions to European options in a regime switching 
economy. Operations Research Letters 33(2005), 581-586. 

Ouhbi, B. Limnios, N., Nonparametric estimation for semi-Markov processes based on its hazard 
rate functions. Stat. Inference Stoch. Process., 2(1999), 151-173. 

Sansom, John and Thomson, Peter, Fitting hidden semi-Markov models to breakpoint rainfall data, 
J. Appl. Probab. 38A (2001), 142157. 

Schweizer, M., A guided tour through quadratic hedging approaches(2001). In E.Jouini, J.Cvitanic, 
& M.Musiela (Eds.), Option pricing interest rates and risk management. Cambridge: Cambridge 
University Press, 538-574. 


12 



